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In our study we consider nonlinear, power-law field-dependent electrical permitivity and magnetic 
permeability and investigate the time-dependent Maxwell equations with the self-similar Ansatz. 
This is a first-order hyperbolic PDE system which can conserve non-continuous initial conditions 
describing electromagnetic shock-waves. Besides shock-waves other interesting solutions (e.g. with 
| localized compact support) can be found with delicate physical properties. Such phenomena may 

happen in complex materials induced by the planned powerful Extreme Light Infrastructure(ELI) 
^ ■ laser pulses. 
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INTRODUCTION 



Wave propagation in non-linear media is a fascinating field in physics with large literaturefjj . To study such 
effects diverse non-linear partial differential equations(PDEs) have to be investigated with various methods. The 
most investigated non- linear wave propagation phenomena are the solitary waves, usually based on the non- linear 
Schrodinger or sine-Gordon or KdV equations. On the other side there are many more, not-so well-known non-linear 
wave equations exist which have delicate properties such as, shock-waves, solutions with continuous compact support 
and so on. Such equations are the various Euler or non-linear heat conduction equations. To investigate if a system 
has such properties one of the most powerful analytical tool is to apply the self-similar Ansatz which may describe 
dispersive solutions with reasonable physical interpretation. The validity of such solutions is very wide in continuum 
mechanics and mostly used to study shock- waves and other fluid dynamical problems 0-0). 

In one of our former studies we investigated the paradox of heat conduction with a new kind of time-dependent 
Cattaneo heat conduction law [j| and found physically reasonable solutions with compact support. As a second 
independent study we presented three dimensional analytical results for the Navier-Stokes equations . The properties 
of the self-similar solution will be analyzed later. 
| It is well-known that from the four Maxwell field equations combining with the two constitutive relations a linear 
■^j- ■ second-order hyperbolical wave equation can be derived for the field variables. In such cases the constitutive equations 
contain only linear relations for the electrical permitivity and for the magnetic permeability. The theory of the 
electromagnetic wave propagation can be found in various studies Q 

Under non-linear Maxwell equation most of the people mean the nonparaxial nonlinear Schrodinger equation 
(NNSE) which is derived from the Helmholtz equation including the Kerr media where the relative dielectric permi- 
tivity is well described e r = n 2 = jiq + 5Nh{E y ). Here no is the linear contribution to the total refractive index n. In 
sufficiently slow media, where the characteristic response time of the nonlinearity is much greater than the temporal 
period of the field oscillations one has 8nl ~ 2rign2 < E 2 > where ri2 is the Kerr coefficient and <> denotes the time 
^ ' averaging over many optical cycles. 

There are large number of studies available where the NNSE is analytically (or numerically) solved and the obtained 
solitons are analyzed. Additional literature can be found in |8|. 

When ultrashort intense laser pulses propagate in a media then there is an intensity dependence of the group 
velocity which leads to the phenomena of self-steepening and optical shock wave formation. It means that the peak 
of the pulse is slowed down more than the edge of the pulse, leading to steepening of the trailing edge of the pulse. 
The envelope becomes steeper and steeper. If the edge becomes infinitely steep, it is said to form an optical shock 
wave. Self-steepening has been described by various authors In non-linear media optical beams can suffer 

self-trapping where the wave equation is solved with the displacement field of D = eoE + rjE 3 [l5|, [l6| . 
Most of the authors consider electromagnetic shock-waves in this sense. 

We however follow a different way, in our recent study we consider nonlinear, power-law field-dependent electrical 
permitivity and for the magnetic permeability and investigate the last two time-dependent Maxwell equations with the 
self-similar Ansatz. This is now a first-order hyperbolical non-linear PDE system which can conserve non-continuous 
initial conditions describing electromagnetic shock-waves. To our knowledge the question of electromagnetic shock- 
wave is only laconically shortly mentioned in a well-known physic textbook (l7j . There is a small number of publi- 
cations exist (most of them are from Russian authors) about electromagnetic shock-wave propagations is anisotropic 
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magnetic materials where the direct Maxwell equations [18| are solved. Unfortunately, not with our method. 

Besides shock-waves other interesting solutions (e.g. with compact support) can be found with delicate physical 
properties. Such effects might happen in complex materials which could be induced by powerful laser pulses which 
will be available in the planned Extreme Light Infrastructure (ELI). 



II. THEORY AND RESULTS 



Let's start with the well-known four Maxwell equations for the fields: 



V-D = p VB = 

VxE = -f VxH^+J (1) 

at at 

where E,B are the electric and magnetic fields D,H are the electric displacement and magnetizing fields p is the 
electric charge and J is the current density which is zero in insulator media. The closing constitutive relations are 
the following: 

D = eE B = //H (2) 

where e is the electrical permitivity and fi is the magnetic permeability. For non-isotropic linear materials e and /i 
are second order tensors. (D a — t a pEp and H a = n a pBp). For the linear and isotropic materials cases these 
are pure real numbers. The most general linear relation for the constitutive equations is however the following 

D a ( x ,t)=J2 I d 3 *' I dt'e a ^,t')E^ -x,i-i')- (3) 

(In Eq. 3 the D a — ¥ B a and Ep — > Hp interchange is still valid.) This equation means non-locality both in space 
and time. The later can be addressed as memory effects, too. The Fourier transform of the electrical permitivity is 
the frequency dependent dielectric function, which attracts much interest. The crucial symmetry properties can be 
expressed via the Kramers-Kronig formula 17], 19] which dchne the relation between the real and the imaginary part. 

Equations of (Q][2]) are enough to derive the well-known second-order liner hyperbolic wave-equation for the field 
variables, which can be found in any electrodynamics textbook 19]. However, this classical calculation is based on a 
numerical trick, an additional radial derivation is done, which also means that the first derivatives of the fields are 
continuous and small. But the original Maxwell equations are of first order both in time and space, therefore the 
initial conditions need not be continuous. 

According to the basic book of Zel'dovich Q to describe the propagation of large disturbances (non-continuous 
tears, shock-waves) in a media the first order hyperbolic Euler and continuity equations have to be applied, or with 
other words these equation also have f(x — at) traveling wave solutions with a velocity of c which is larger then 
the propagation of sound. The speed of sound however enters the gasdynamic equations as the velocity of small 
propagation of small disturbances hence the usual second-order wave equation. In this language we may speak about 
two different kind of wave equations or wave propagation phenomena. The Maxwell equations should be considered 
for large disturbances and the second order wave equation for the small, sinusoidal disturbances. This is the basic 
idea in our recent study, and we want to study the existence of shock waves. We follow this analogy and apply 
non-linear material laws and solve directly the first-order hyperbolic Maxwell equations for propagation and instead 
of the second-order wave equation. 

Maxwell equations in vacuum are linear in the fields of B and E. Many hundred telephone conversation can 
propagate parallel on a single microwave link without any distortion. Another experimental evidence of the linearity 
is the idea of linear superposition. In optics white light is refracted by a prism into the color of the rainbow and 
recombined into white light again. There are of course, circumstances when non-linear effects occur in magnetic 
materials or in crystals responding to intense laser beams like frequency doubling. 

Our non-linear Maxwell equation is, however defined in a complete different way, namely through the following 
non-linear (material) or constitutive equations 

/i(H) = aH<? e(E) = bW (4) 



where all the four free parameters (a,b,q,r) are real numbers (for physical reasons e(E) • /i(H) > ) and the a 
and b are present to fix the physical dimensions. (Such power law dependence of material constants are well known 
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in different flow problems like in heat propagation [3J where the heat conduction constant can have temperature 
dependence n ~ T v .) Note, that via this relation we define space and time dependent (dynamical) material equations 
which are still local in space and time. 

We know from special relativity that the speed of light in vacuum is the the largest available wave propagation which 
can carry physical information, and can be evaluated from the electric and magnetic properties as well c 2 = 1/ (/io£o)- 
(The zero subscript stands for vacuum.) Permeability and permitivity are not fully independent from each other. 
This formula is slightly modified for any additional media like (? m — l/(/io/^m£oem) where the subscript m stands 
for the media. It is also clear that any stable electromagnetic wave propagation speed in media is always less than 
the speed of light in vacuum (c m < c) but for a short time quick particles (usually charged) can propagate quicker 
than the local speed of light producing Cherenkov radiation. Therefore in our calculations we will use the following 
relation: 

v(H) = am e (#) = _j_. (5) 

note that now the propagation speed of the electromagnetic signal has an upper bound which is c. (From now on 
we will consider one spatial coordinate and neglect the vectorial notation.) Now we reduced the number of free 
parameters to two. In electromagnetic wave propagation the role of e and /x are symmetric, however we use this 
relation because of the existence of J in the last Maxwell equation. We will see later on that with this choice the 
ordinary differential equation which is obtained from the third Maxwell equation can be integrated and the solutions 
became more transparent. 

For the current density we apply the differential Ohm's law: 

J = <tE (6) 

where a is the conductance of the media - and can be a second rank tensor in crystals or a highly non-linear 
field dependent quality like the permeability or the susceptibility a = hEi p . In a transition- metal oxide it can be a 
a fa (l/E)sinh(E) function go|. 

(Our physical intuition says that only some integers (=f1,=f2) and some rational numbers (^1/2,^2/3) will be 
crucially interesting.) 

For the sake of simplicity we consider the following one dimensional wave propagation problem 

E=(0,E v (x,t),0) H = (0,0, H z (x,t)) (7) 

which means linearly polarized electric filed in y direction with x coordinate dependence and linearly polarized 
magnetic field in z direction with x coordinate dependence only. Now the last two Maxwell equations are: 



dE y _ dB z dH z _ dD y 

From basic textbooks the form of the one-dimensional self-similar Ansatz is well-known 

T(x,t)=t-<*f(^);=t- a f(r ] ). (9) 

where T(x,t) can be an arbitrary variable of a partial differential cquation(PDE) and t means time and x means 
spatial dependence. The similarity exponents a and j3 are of primary physical importance since a represents the rate 
of decay of the magnitude T(x, t), while /3 is the rate of spread (or contraction if /3 < ) of the space distribution as 
time goes on. Solutions with integer exponents are called self-similar solutions of the first kind (and sometimes can 
be obtained from dimensional analysis of the problem) . The above given Ansatz can be generalized considering real 
and continuous functions a(t) and b(t) instead of t a and t@. The applicability of self-similarity also means that no 
characteristic time scale (like a relaxation time) is presented in the investigated system. 

The most powerful result of this Ansatz is the fundamental or Gaussian solution of the Fourier heat conduction 
equation (or for Fick's diffusion equation) with a = j3 = 1/2. Applicability of this Ansatz is quite wide and comes up 
in various transport systems 

For our problem we consider the following Ansatzc: 

E y (x,t) = t- a f (J) := t- a f(v); H z (x,t) = r s g := t^g^). (10) 

Where a, /3, 5 are three real number which are (at this point of the model) independent from each other. The functions 
/(?]) and g(r)) are the shape functions of the problem. 
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Combining Eq. (fT0|) together with Eq. ([5]), (HJ and inserting into the original last two Maxwell equations we get 
the following system: 



Having done the derivations we arrive at the next ordinary differential equation(ODE) system 

/' = a(q+l)[Sg^ +1 -g^ V 0\ 
-9 = -^l(q + l)9 q f + q(q + l)9 q - 1 g / fv+(q + V9 q fv} 



(11) 



(12) 



where prime means derivation with respect to rj. Note, that if the following universality relations are held (5=1 and 
j3 = —q — 1) the first equation is a total derivative and can be integrated resulting 



f = a(q+l) V g«+ 1 . 



(13) 



This fixes the connection between the electric and magnetic fields. From the second ODE we get that a = — 1 and 
p = 1 should be. This means that our media should not have any conductivity for self-similar solutions. Inserting Eq. 
(|13jl into the second equation of (|12|) we arrive to our final expression: 
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2{q + l) 2 V g 2 i +1 + (g + 1) 3 7? 2 g 2 g + h 
l + q(q+ l) 2 f] 2 g 2 i 



(14) 



where light velocity c is fixed to unity remaining q and h the final two free parameters. The role of q is more 
important than the role of h. Note, that now different real q values mean different exponents for magnetic permeability 
representing different physical material properties and different physics. 

For general q there is no closed form for the solutions therefore we investigated the direction field of the ODE 
which helps us to gain global information for the solutions. Note that Eq. (|14p is non-autonomous (depending on rj) 
therefore there is no general theorem for the direction field available. 

It is clear that for q < the denominator can be zero, therefore a singularity can appear where the first derivative 
of f(rf) can be infinite, which means vertical direction field. If a solution with an initial condition meets this field line 
than it stops, and cannot be continued. Note, that if q < —1 the denominator is zero for a g(ri) ~ rf where r < 1, 
which means that the points where the first derivative is singular lie on a concave function. Our numerical experience 
shows that for q < —5/4 all the solutions with an initial condition of /(0) = r where r is a positive real number have 
a compact support with a finite 770 and a final jump from f(r)o) down to zero. 

Figure la presents the direction field for q = —2. Solutions with compact support are clearly presented. This means 
an E y = t~ 1 g(x ■ t) dependence. The compact support of the ODE solution of Eq. (TT4")) means that the solution of 
the original PDE system for E y (x,t) is also compact via the rj = x ■ t in real time and space. The constrain (Eq. 
(fl3|) ) between f and g also dictates the same compact support for the H z (x,t) field as well. Outside these time and 
x coordinate ranges we may fix the values of E y (x,t) and H z (x,t) identically to zero which are also solutions of the 
last two Maxwell equations. In this way we got shock-wave solutions for the above mentioned trajectories. 

If q > —1 than all the solutions are continuous. For — 1 < q < 1 all the initial conditions /(0) = r where r is a 
positive real number give blow-up solutions which means compact domain and infinite range. For q > 1/2 all the 
initial conditions f(p) — r where p and r are positive real numbers give solutions which have no limit in domain or in 
range. Fig. lb presents such solutions for q = 2. 

There are two analytic solutions available: The q = — 1 limit can be taken and gives the g' = ODE which has the 
constant as a trivial solution. 

The equation has a second explicit closed solution for q=0 which is the following: 
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2e r ' 



\pKErr(r\) + —h\fiiErr{r\) + < 



(15) 



here Err(rf) means the usual error function [21| . Note that q=0 gives back the usual linear wave equation which 
might have such a continuous solution which goes to infinity for large 77. So this is not much new. 

Asymptotic analysis of the denominator of Eq. (|14|) can be be considered too. Unfortunately, for 1 >> q(q+l) 2 i] 2 g 2q 
there is no closed solution for the remaining ODE for general q. There is an analytic solution for q = 1/2 only 
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FIG. 1: The direction field of Eq. (|14[) a) for g = —2 and b) g = +2. Different lines shows solutions for different initial 
conditions. 



= (3(-r/ 3 ) 2 / 3 e 9 / 8r ' 3 )/[23 2/ Vr(2/3) - 23 2/ Vr(2/3, -9/8ry 3 ) - 3ci(-?7 3 ) 2 / 3 ] where T means the gamma function 
[21j . This function has a singularity between 1/2 and 1. 

For the other case when y(y + l) 2 n 2 g 2q >> 1 there the remaining ODE is 5' = (779 + 2y + rj)/{r)q) with the 
analytic solution of g = (rj 2 ^ q ci +f?(g'+l))/(9— 2) with arbitrary q. The integration constant is notated with ci as usual. 

To find physically reasonable solutions we calculate the Poynting vector which gives us the energy flux (in W/m 2 ) 
of an electromagnetic field. Unfortunately, there are two controversial form of the Poynting vector in material based 
on the Abraham or the Minkowski formalism, a detailed description can be found in [22[. We use the classical form 
of Poynting vector 

S = E x H = t- a - s fg = r a - s a{q + l)rjg q+2 . (16) 

Note that for q < —2 the f° u Sdrj is finite which is a good result. The spacial integral of the pointing vector 

~ J™* w9 q+2 (w*)dx is finite for all time. However the time integral of ~ J^"* ■^■g q+2 (-^-)dt can be problematic for 
small t, and depends on the concrete form of g. 

Another method to classify if the solutions are physical would be the total energy of the fields in a finite volume. 
For linear electrodynamics the energy density is defined as follows W = |(E • D + B • H) However, even this formula 
is not unique there are various non-linear electrodynamics by Born 23J or by Rafelski [2~ij based on the Lagrangian 
density which contains additional terms. Kotel'nikov [25| generalized the Born model and suggested an infinite series 
of Lorentz and Poincare-invariant nonlinear versions of the Maxwell equations. 

Our approximation to describe the permeability and permittivity Eq. ([5]) is just one way to a non-linear model. 
Another physically tenable description for the constitutive equations could be a series expansion like e(E) = 1 + 
aE + bE 2 ...) where the linear term is responsible for the so-called Pockels or electro-optical effect and the quadratic 
term is the for the Kerr effect (a and b are constants to fix the proper dimension). For optically important materials 
fi = const is the right choice, non-linear magnetic properties play a significant role only for plasmas. Unfortunately, 
our well-established Ansatz Eq. (jlOl) does not applicable directly to such power series. 
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III. SUMMARY 



We introduced a power law magnetic field dependence of the magnetic permeability and investigated the corre- 
sponding non-linear Maxwell filed equations with the self-similar Ansatz. If the exponent of the field dependence of 
the permeability is smaller that minus one than compact, shock- wave like solutions are obtained which might have 
some importance in laser matter interactions. The work was supported by the Hungarian HELIOS project. 
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